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Abstract 

The aim of this note is to provide a pedagogical survey of the recent works [30, 31] 
concerning the local behavior of the eigenvalues of large complex correlated Wishart ma¬ 
trices at the edges and cusp points of the spectrum: Under quite general conditions, the 
eigenvalues fluctuations at a soft edge of the limiting spectrum, at the hard edge when 
it is present, or at a cusp point, are respectively described by mean of the Airy kernel, 
the Bessel kernel, or the Pearcey kernel. Moreover, the eigenvalues fluctuations at several 
soft edges are asymptotically independent. In particular, the asymptotic fluctuations of 
the matrix condition number can be described. Finally, the next order term of the hard 
edge asymptotics is provided. 


1 The matrix model and assumptions 

Consider the N x N random matrix defined as 

= AxjvXArXy (1.1) 

where Xjv is an A x n matrix with independent and identically distributed (i.i.d.) entries 
with zero mean and unit variance, and Xw is a n x n deterministic positive definite Hermitian 
matrix. The random matrix Mw has N non-negative eigenvalues, but which may be of 
different nature. Indeed, the smallest N — min(n, N) eigenvalues are deterministic and all 
equal to zero, whereas the other min(n, N) eigenvalues are random. The problem is then to 
describe the asymptotic behavior of the random eigenvalues of Mw, as both dimensions of 
Xw grow to infinity at the same rate. Let us mention that the n x n random covariance 
matrix 

Mw = ^xj^/'xyxwx;^/', 

which is also under consideration, has exactly the same random eigenvalues as Mw, and hence 
results on the random eigenvalues can be carried out from one model to the other immediately. 
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The global behavior of the spectral distribution of M^r in the large dimensional regime 
is known since the work of Marcenko and Pastur [42], where it is shown that this spectral 
distribution converges to a deterministic probability measure /i that can be identified. In this 
paper, we will be interested in the local behavior of the eigenvalues of M^r near the edge 
points and near the so-cjdled cusp points of the support of fj,. The former will be called the 
extremal eigenvalues of M^v • 

The random matrices M^v and M^r are ubiquitous in multivariate statistics [4], math¬ 
ematical finance [39, 46], electrical engineering and signal processing [18], etc. Indeed, in 
multivariate statistics, the performance study of the Principal Component Analysis algo¬ 
rithms [36] requires the knowledge of the fluctuations of the extremal eigenvalues of M^v. 
In mathematical finance, M^r represents the empirical covariance matrix obtained from a 
sequence of asset returns. In signal processing, M^r often stands for the empirical covariance 
matrix of a spatially correlated signal received by an array of antennas, and source detec¬ 
tion [12, 38] or subspace separation [45] algorithms also rely on the statistical study of these 
extremal eigenvalues. 

In this article, except when stated otherwise, we restrict ourselves to the case of complex 
Wishart matrices. Namely, we make the following assumption. 

Assumption 1. The entries of are i.i.d. standard complex Gaussian random variables. 

Concerning the asymptotic regime of interest, we consider here the large random matrix 
regime, where the number of rows and columns of M^v both grow to infinity at the same pace. 
More precisely, we assume n = n{N) and n, A —)> oo in such a way that 

Tl 

lim ^ = 7 e (0,oo) • (1-2) 

N^oo iV 

This regime will be simply referred to as A —)• oo in the sequel. 

Turning to Xtv, let us denote by 0 < Ai ^ ^ An the eigenvalues of this matrix and let 

1 ” 

= (1-3) 

J=1 

be its spectral measure. Then we make the following assumption. 

Assumption 2. 

1. The measure weakly converges towards a limiting probability measure as A —?• oo, 
namely 

n 

i=i 

for every bounded and continuous function /. 

2. For A large enough, the eigenvalues of stay in a compact subset of (0,-|-oo) inde¬ 
pendent of A, i.e. 

0 < liminfAi, sup An < -|-oo. (1-4) 

N^oo jy 


2 



Under these assumptions, a comprehensive description of the large N behavior of the 
eigenvalues of M^r can be made. To start with, we recall in Section 2 some classical results 
describing the global asymptotic behavior of these eigenvalues, as a necessary step for studying 
their local behavior. We review the results of Marcenko-Pastur [42] and those of Silverstein- 
Choi [50], which show among other things that the spectral measure of Mjv converges to a 
limit probability measure //, that /x has a density away from zero, that the support of /x can be 
delineated, and that the behavior of the density of /x near the positive endpoints (soft edges) 
of this support can be characterized. We moreover complete the picture by describing the 
behavior of the limiting density near the origin when it is positive there (hard edge), and also 
when it vanishes in the interior of the support (cusp point). The latter results are extracted 
from [31]. 

Next, in Section 3 we turn to the eigenvalues local behavior. More precisely, we investigate 
the behavior of the random eigenvalues after zooming around several points of interest in the 
support, namely the soft edges, the hard edge when existing, and the cusp points. In a word, 
it is shown in the works [30, 31] that the Airy kernel, the Bessel kernel, and the Pearcey 
kernel describe the local statistics around the soft edges, the hard edge, and the cusp points 
respectively, provided that a regularity condition holds true. In particular, the extremal 
eigenvalues fluctuate according to Tracy-Widom laws. 

In Section 4, we provide sketches of proofs. We first recall an important expression for 
the kernel K^v associated to the (random) eigenvalues of and then outline how one can 
prove asymptotic convergence towards Airy, Pearcey or Bessel kernels by zooming around the 
points of interest: either a soft edge, a cusp point or the hard edge. 

In Section 5, we provide a list of open questions, directly related to the results of the 
paper. 

Acknowledgements. WH is pleased to thank the organizers of the Journees MAS 2014 
where the project of this note was initiated. During this work, AH was supported by the 
grant KAW 2010.0063 from the Knut and Alice Wallenberg Foundation. The work of WH 
and JN was partially supported by the program “modeles numeriques” of the French Agence 
Nationale de la Recherche under the grant ANR-I2-MONU-0003 (project DIONISOS). 

2 Global behavior 

Since the seminal work of Marcenko and Pastur [42], it is known that under Assumptions I 
and 2 the spectral measure of Mw almost surely (a.s.) converges weakly towards a limiting 
probability measure /x with a compact support. Namely we have 

^Tr/(Mjv) j /(x)/x(dx) (2.1) 

for every bounded and continuous function /. As a probability measure, /x can be character¬ 
ized by its Cauchy transform: this is the holomorphic function defined by 

m{z) = f -- /x(dA), z G C+ = {z G C : Im( 2 ;) > O}, 

J z — X 
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and which takes its values in C_ = {z € C : Ini( 2 ;) < 0}. More precisely, for any open 
interval / C M with neither endpoints on an atom of /x, we have the inversion formula 

u,(I) = -lim [ lm(m(x + ie))dx. 

TT e^O Jj^ ’ 

For every z G C+, the Cauchy transform m{z) of /x happens to be the unique solution m G C_ 
of the fixed-point equation 

where 7 and v were introduced in ( 1 . 2 ) and Assumption 2 . 

Moreover, using the free probability terminology [2, 32], the limiting distribution /x is also 
known to be the free multiplicative convolution of the Marcenko-Pastur law (2.3)-(2.4) with 
V, and equation ( 2 . 2 ) is a consequence from the subordination property of the multiplicative 
free convolution [19]. 

For example, in the case where v = 81 , which happens e.g. when Sjv = In, this equation 
has an explicit solution and the measure ^ can be recovered explicitly. 

/x(dx) = (1 - 7 )+(5o +/c(a^)dx, (2.3) 

where x^ = max(x, 0 ) and the density p has the expression 

p{x) = {h - x){x - a) l[„,t,](x), a =(1-^)2, h = {l + ^f. (2.4) 

This is the celebrated Marcenko-Pastur law. 

When V has a more complicated form, it is in general impossible to obtain an explicit 
expression for p,, except in a few particular cases. Nonetheless, it is possible to make a 
detailed analysis of the properties of this measure, and this analysis was done by Silverstein 
and Choi in [50]. These authors started by showing that = fnix) exists for 

every x G M* = M — {0}. Consequently, the function m{z) can be continuously extended to 
C+ U M*, and furthermore, /x has a density on M* defined as p{x) = — 7 r“^Im (m(x)). We still 
have the representation 

p{dx) = {1 - 80 + pix)dx (2.5) 

with this new p, making p{x)dx the limiting distribution of the random eigenvalues of Mtv. 
As is common in random matrix theory, we shall refer to the support of p{x)dx as the bulk; 
we will denote (with a slight abuse of notation) its support by Supp(/?). Silverstein and Choi 
also showed that p is real analytic wherever it is positive, and they moreover characterized 
the compact support Supp{p) following the ideas of [42], More specifically, one can see that 
the function m{z) has an explicit inverse (for the composition law) on m(C+) defined by 

g{m) = — + 7 / ^ . i^(dA), (2.6) 

m 7 1 — m\ 

and that this inverse extends to C_ UT> and is real analytic on D, where D is the open subset 
of the real line 

D = {x G M : X / 0, x~^ ^ Supp(zx)}. (2.7) 
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Figure 1: Plot of : L) —M for 7 = 0.1 and v = 0.75i + 0.3(53. In this case, 
D = (—00,0) U (0, |) U (|,1) U (1,00). The two thick segments on the vertical 
axis represent Supp(/?). 


It was proved in [50] that 

M — Supp(p) = {g{m) : m & D, g'{m) < O}. 


An illustration of these results is provided by Figures 1 and 2. 

Of interest in this paper are the left edges, the right edges and the cusp points of Supp(/9). 
A left edge is a real number a satisfying for every (5 > 0 small enough 

a+5 


ra r<x+d 

/ p{x)dx = 0, / p{x)(ix > 0 . 

' a—(5 J a 


A right edge is a real number a satisfying for every 5 > 0 small enough 

ra ra+S 


' a—(5 


p{x)dx > 0, 


p{x)dx = 0 . 


A cusp point is a real number a such that p{a) = 0 and, for every <5 > 0 small enough, 

pa paH-5 

/ p{x)dx > 0 and / p{x)dx > 0 . 

J a—5 J a 

Of course all edges and cusp points are positive numbers, except perhaps the leftmost 
edge. When the leftmost edge is the origin, it is common in random matrix theory to refer 
to it as the hard edge. In contrast, any positive edge is also called a soft edge. 

The results of [50] summarized above show that the study of the map g on the closure 
D oi D provides a complete description for the edges and the cusp points. First, a right 
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Figure 2: Plot of the density p in the framework of Figure 1. 


edge is either a local minimum of g reached in D, or belongs to g{dD), which means there is 
c G dD = D\D such that 9{x) exists, is finite, and equals to that edge. In the 

former case, p{x) behaves like a square root near the edge. 

Proposition 2.1. If a is a right edge, then either there is a unique c € D such that 

g{c) = a, g\c) = 0, ff"(c) > 0, (2.8) 

or a £ g{dD). In the former case, we have 



Conversely , if c £ D satisfies (2.8), then a is a right edge and (2.9) holds true. 

The case where an edge lies in g{dD) turns out to be quite delicate. In the forthcoming 
description of the eigenvalues local behavior near the edges, we shall restrict ourselves to the 
edges arising as local minima of g, see also Section 5 for further discussion. Notice also that 
if is a discrete measure, as exemplified by Figures 1 and 2, then g is inhnite on dD and in 
particular a right edge cannot belong to g{dD): the right edges are in this case in a one-to-one 
correspondence with the local minima of 5 on H. 

The situation is similar for the soft left edges, except that they correspond to local maxima. 

Proposition 2.2. If a > 0 is a left edge, then either there is a unique c £ D such that 

g{c) = a, g'{c) = 0, g"(c) < 0, (2.10) 
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or a G g{dD). In the former ease, we have 

1 / O \ 1/^ 

Pi^) = - + 0(1)) , x-^a+. (2.11) 

Conversely, if c ^ D satisfies (2.10), then a is a right edge and (2.11) holds true. 

Propositions 2.1 and 2.2 have been established in [50]. We state below their counterparts 
for the hard edge and a cusp point. 

The hard edge setting turns out to be similar to the soft left edge one, except that c is 
now located at infinity, and p{x) behaves like an inverse square root near the hard edge. More 
precisely, observe that the map g is holomorphic at oo and g{oo) = 0, in the sense that the 
map z I—)• g(l/z) is holomorphic at zero and vanishes at z = 0. We also denote by g'{oo) and 
g''{oo) the first and second derivatives of the latter map evaluated at z = 0. 

Proposition 2.3. The bulk presents a hard edge if and only if 

g{oo) = 0, g'{oo) = 0, /(oo) < 0, (2.12) 

or equivalently i/ 7 = 1. In this ease, we have 

p{x) = -( —(1 + 0(1)) , x^0+. (2.13) 

vr \-g”{oo)J 

More precisely, we have the explicit formulas g'{oo) = 1 —7 and g”{oo) = —27 f A~^n(dX). 
In particular the statement 5''(00) < 0 is always true, and so is g(oo) = 0 as explained above; 
we included them in (2.12) to stress the analogy with (2.10). 

A simple illustration of Propositions 2.1 to 2.3 is provided by the Marcenko-Pastur law. 
From (2.4), one immediately sees that p(x) rsj (t’ — x)^/^ as X —)• b_, and that a similar square 
root behavior near a holds if and only if a > 0, that is 7 7^ 1. If 7 = 1, i.e., a = 0, then 
p(x) ~ x“^/^ as X —)• 0+ instead. 

We now turn to the cusp points. Those who will be of interest here correspond to inflexion 
points of g where this function is non decreasing. Moreover, a cubic root behavior for the 
density p{x) is observed near such a cusp point, hence justifying the terminology (we recall 
cusp usually refers to the curve defined by = x^). 

Proposition 2.4. Let a be a cusp point, set c = m{a) and assume c £ D. Then 

5(c) = a, g'{c) = 0, /(c) = 0, and /'(c) > 0. (2.14) 

Moreover, 

^ (/^) k “ + ’ x^a. (2.15) 

Conversely, if c £ D satisfies y'(c) = y"(c) = 0, then the real number a = y(c) is a cusp point, 
/'(c) > 0 and (2.15) holds true. 

Propositions 2.3 and 2.4 appear in [31]. Proposition 2.4 is illustrated in Figures 3 and 4. 
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Figure 3: Plot of g for 7 ~ 0.336 and v = 0.7(5i + 0.3(53. The thick segment on the vertical 
axis represents Supp(/r). The point a is a cusp point. 

3 Local behavior 

The study of the eigenvalues local behavior of random matrices is a central topic in random 
matrix theory. When dealing with large Hermitian random matrices, it is recognized that the 
local correlation of the eigenvalues around an edge where the density vanishes like a square 
root should be described by a particular point process involving the Airy kernel (see below), 
whose maximal particle’s distribution is known as the Tracy-Widom law. For instance, this 
has been established for unitary invariant random matrices and for Wigner matrices as well, 
see e.g. the surveys [21, 27] and references therein. Similarly, the Bessel kernel is expected to 
describe the fluctuations around an hard edge where the density vanishes like an inverse square 
root, and the Pearcey kernel around a cusp with cubic root behavior. Let us also mention 
that the sine kernel is expected around a point where the density is positive, and that more 
sophisticated behaviors have also been observed in matrix models where the density vanishes 
like a rational power of different order, but we will not further investigate these aspects here. 
Another interesting feature not covered by this survey is the study of the random eigenvectors, 
see e.g. [10, 13]. 

The purpose of this section is to present the central results of [30, 31], where such typical 
local behaviors arise for the complex correlated Wishart matrices under consideration at every 
edges and cusp points satisfying a certain regularity condition. 

In fact, the free parameter family {i^n) may have a deep impact on the limiting local fluc¬ 
tuations and one may not recover the expected fluctuations without further conditions. A hrst 
manifestation of this phenomenon is the Baik-Ben Arous-Peche (BBP) phase transition, that 
we present in Section 3.1. In a nutshell, this phase transition yields that slight variations on 
the family (vn) may modify the fluctuations at a soft edge and may no longer be described by 
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Figure 4: Plot of the density of n in the framework of Figure 3. 


the Tracy-Widom law. Such phenomenas motivate the introduction of a regularity condition 
which essentially rules out this kind of behaviors. 

In Section 3.2, we provide the existence of finite N approximations of the edges or cusp 
points under study and satisfying the regularity condition; the reader not interested in these 
precise definitions may skip this section. 

Next, in Section 3.3 we introduce the Airy kernel, the Tracy-Widom law, and state our 
results concerning the soft edges. In Section 3.4 , we introduce the Bessel kernel and describe 
the fluctuations at the hard edge. As an application, we provide in Section 3.5 a precise 
description for the asymptotic behavior of the condition number of Mjv- Finally, in Section 
3.6 we introduce the Pearcey kernel, and state our result concerning the asymptotic behavior 
near a cusp point. 

3.1 The BBP phase transition and the regularity assumption 

First, assume Sat is the identity matrix, so that the (limiting) spectral distribution of S^r 
is (5i, and hence the limiting density p{x) is provided by (2.4). If Xmax stands for the maximal 
eigenvalue of Mat, then it has been established that Xmax converges a.s. towards the right 
edge and fluctuates at the scale according to the Tracy-Widom law [34]. Next, following 
Baik, Ben Arous and Peche [7], assume instead Sat is a finite rank additive perturbation of the 
identity, meaning that the rank of the perturbation is independent on N. Thus we still have 
V = 5i and the limiting density p{x) remains unchanged. They established that if the strength 
of the perturbation is limited, then the behavior for Xmax is the same as in the non-perturbed 
case, see [7, Theorem 1.1(a), k = 0]. On the contrary, if the perturbation is strong enough, 
then Xmax converges a.s. outside of the bulk and the fluctuations are of different nature, see [7, 
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Theorem 1.1(b)]. But in this case, one can consider instead the largest eigenvalue that actually 
converges to the right edge and show that the Tracy-Widom fluctuations still occur (this is 
a consequence of Theorem 1 below). However, they also established there is an intermediary 
regime, where Xmax converges a.s. to the right edge and fluctuates at the scale but 

not according to the Tracy-Widom law, see [7, Theorem 1.1(a), k > 1], hence leaving the 
random matrix universality class since the right edge exhibits a square root behavior. Here 
the fluctuations are actually described by a deformation of the Tracy-Widom law, but in the 
general XItv’s setting much exotic behaviors must be expected. 

In conclusion, although the eigenvalues global behavior only depends on the limiting 
parameters v and 7, the local behavior is quite sensitive in addition to the mode of convergence 
of the spectral measure of Sat to its limit u. In order to obtain universal fluctuations in 
the more general setting under investigation, it is thus necessary to add an extra condition 
for the edges, and actually for the cusp points too. A more precise consideration of the non- 
universal intermediary regime considered by Baik, Ben Arous and Peche reveals that, if we 
write the right edge as g{c), see Proposition 2.1 and the comments below, then some of the 
inverse eigenvalues of 51 at converge towards c. Recalling the Aj’s stand for the eigenvalues of 
Hat, this motivates us to introduce the following condition. 


Definition 3.1. A real number c satisfies the regularity condition if 


n 

lim inf min 

N^OO j = l 



(3.1) 


Moreover, if c satisfies the regularity condition, we then say that g{c) is regular. 

Remark 3.2. Propositions 2.1 and 2.2 tell us that every soft edge reads g{c) for some c £ D. 
In fact, since 5(0) = -|-oo and Supp(/x) is compact, necessarily c 7^ 0. If we moreover assume 
the soft edge to be regular then, since by definition D = {x G M : x 7^ 0, x~^ 0 Supp(z^)} and 
because z^a^ converges weakly to i', necessarily c £ D. In particular. Propositions 2.1 and 2.2 
yield that at a regular soft edge the density show a square root behavior. As regards the hard 
edge, the analogue of the regularity condition turns out to be lim inf at Ai > 0 and is therefore 
contained in Assumption 2. 


Remark 3.3. We show in [30] that, if 7 > 1, then the leftmost edge a is always regular. 
Namely there exists c £ D which is regular such that a = g{c). In fact, we have c < 0. 

Before we state our results on the eigenvalues local behavior around the regular edges or 
the cusp points, we now provide the existence of the appropriate scaling parameters we shall 
use in the later statements. 


3.2 Consequences of the regularity condition and finite N approximations 
for the edges and the cusp points 

Recall from Section 2 that the Cauchy transform of the limiting eigenvalue distribution g 
of Mat is defined as the unique solution m £ C_ of the fixed-point equation (2.2). We 
now consider the probability measure induced after replacing (7, z^) by its finite horizon 
analogue in this equation (we recall z^az was introduced in (1.3)). Namely, let 

be the probability measure whose Cauchy transform is defined as the unique solution m £ C_ 
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of the fixed-point equation 

The probability measure should be thought of as a deterministic approximation of the 
distribution of the eigenvalues of Mjv at finite N, and is referred to as the deterministic 
equivalent of the spectral measure of Mv- The measure /iv reads 

( Tl \ ~^ 

1 - ^ j 5o-h/9Ar(a:)dx, 

and one can apply all the results stated in Section 2 to describe pN, after replacement of g 
with 

9n{z) = 1 + ^ j I 

Recalling D has been introduced in (2.7), the following proposition encodes the essential 
consequence of the regularity condition. 

Proposition 3.4. If c £ D satisfies the regularity condition (3.1), then there exists <5 > 0 
such that gN is holomorphic on {z £ C : \z — c\ < 6} C D for every N large enough and 
converges uniformly towards g there. 

It is an easy consequence of Montel’s theorem. Now, if a sequence of holomorphic function 

hjsf converges uniformly to a (holomorphic) function h on an open disc, then a standard result 

(k) 

from complex analysis provides that the A:th order derivative hj^' also converges uniformly 
to h^^^ there, for every k 1. Moreover, Hurwitz’s theorem states that, if h has a zero c 
of multiplicity I in this disc, then hjsi has exactly I zeros, including multiplicity, converging 
towards c as —)• oo. 

Thus, as a consequence of the previous proposition, by applying Hurwitz’s theorem to g'j^ 
(and the symmetry g'j^{z) = g'j^^{z)), it is easy to obtain the following statement. 

Proposition 3.5. Assume c £ D satisfies the regularity condition (3.1) and moreover 

9'(0 = 0, g"{c) < 0, resp. g"{c) > 0. 

Then there exists a sequence {cn), unique up to a finite number of terms, converging to c and 
such that, for every N large enough, we have £ D and 

lim 57 v(cv) = g{c), 5Ar(cAr) = 0 , £iv(cAr) < 0 , resp. g'fi{<^N) > 0 . 

N^oo 

Having in mind Propositions 2.2 and 2.1, this proposition thus states that if one considers 
a regular left (resp. right) soft edge a, and thus a = g{c) with c G H by Remark 3.2, then 
there exists a sequence, unique up to a finite number of terms, of left (resp. right) soft edges 
flv = fl' 7 v(civ) for the deterministic equivalent pj\f converging towards o. These soft edges (ajv) 
are finite N approximations of the edge a, while the cv’s are finite N approximations of the 
preimage c. 

When dealing with regular cusp points, the situation is slightly more delicate. The reason 
for this is that if a = g{c) is a regular cusp point, then c is now a zero of multiplicity two 
for gl By applying Hurwitz’s theorem to g'j.^ as above, one would obtain two sequences of 
non-necessarily real zeros for g'j.^ converging towards c. It is actually more convenient to apply 
Hurwitz’s theorem to g'fi instead, in order to get the following statement. 
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Proposition 3.6. Assume c £ D satisfies the regularity condition (3.1) and moreover 

g{c) = 0, = 0, (hence g"'{t) > 0 by Prop. 2.4). 

Then there exists a sequence (cat), unique up to a finite number of terms, converging to c and 
such that, for every N large enough, we have G D and 

lim gN{cN) = 5(c), lim 5iv(civ) = 0, g'l^icN) = 0, 5iv(civ) > 0. 

N^oo N^oo 

Notice that Proposition 3.6 doesn’t guarantee that g'j^{ci\f) = 0. Hence, a cusp point is 
not necessarily the limit of cusp points of the deterministic equivalents /tat. As we shall see in 
Section 3.6, the speed at which goes to zero will actually influence the local behavior 

around the cusp. 

Definition 3.7. Given a soft left edge, resp. right edge, resp. cusp point a which is regular, 
and thus a = g{c) with g'{c) = 0 and g"{c) < 0, resp. g"{c) > 0, resp. g"{c) = 0 and g"'{c) > 0, 
the sequence associated with a is the sequence (cat) provided by Propositions 3.5 and 3.6. 

Equipped with Propositions 3.5 and 3.6, we are now in position to state the results con¬ 
cerning the local asymptotics. 


3.3 The Airy kernel and Tracy-Widom fluctuations at a soft edge 

Given a function K(x, y) from M x M to M satisfying appropriate conditions, one can consider 
its associated determinantal point process, which is a simple point process on M having as 
correlation functions the determinants det[K{yi,yj)]. More precisely, it is a probability dis¬ 
tribution P over the configurations (y*) of real numbers (the particles), namely over discrete 
subsets of M which are locally finite, characterized in the following way: For every k 1 and 
any test function <1> : —)• M, 


E 


Vii 7 ^ ■■■ 7 ^ 


1 k 


^(yi, • ■ ■ ,5fc)det K{yi,yj) dyi---dyk, 


J i,j=l 


where the sum runs over the /c-tuples of pairwise distinct particles of the configuration (yi). 
Hence the correlation between the particles yfs is completely encoded by the kernel K{x,y). 
In particular, the inclusion-exclusion principle yields a closed formula for the gap probabilities 
in terms of Fredholm determinants. Namely, for any interval J C M, the probability that no 
particle lies in J reads 


P 




det 


1 k 


^iyi,yj) , . dyi • • - dyfc , 


J i,j=l 


and the latter is the Fredholm determinant det(/ — K)a 2 (j) of the integral operator acting on 
L^( J) with kernel K(x, y), provided it makes sense. We refer to [33, 35] for further information 
on determinantal point processes. 

Gonsider the Airy point process PAiry which is defined as the determinantal point 
process on M associated with the Airy kernel 


= Ai(x)Ai'(„) - Ai(rtAi'(i) 

X — y 


(3.3) 
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where the Airy function 


1 


OO /^3 


Ai(x) = — / cos 

TT Jo 


+ ux\ du 


is a solution of the differential equation f''{x) = xf{x). 

The configurations (y*) generated by the Airy point process a.s. involve an infinite number 
of particles but have a largest particle ymax- The distribution of ymax is the Tracy-Widom 
law (see e.g. [35, Section 2.2]), and its distribution function reads, for every s G M, 

TAiry(ymax ^ s) = ^Airy ((i/i) H (s, +Oo) = 0^ = det(/ - KAiry)L2(^ ^o) • (3-4) 

Tracy and Widom [51] established the famous representation 

TAiry (?/max ^ s) = exp (^- - s)q{xfdx'^ , 

where q is the Hastings-McLeod solution of the Painleve II equation, namely the unique so¬ 
lution of f”{x) = f{x)^ + xf{x) with boundary condition f{x) ~ Ai(x) as a: —)■ -|-oo. 


Recalling that qm has been introduced in (3.2), we are now in position to describe the 
eigenvalues local behavior around regular soft edges. In the three upcoming theorems, we 
denote by xi ^ ^ the ordered eigenvalues of Ma- We also use the notational convention 

xo = 0 and Xn+i = -|-oo. 

Theorem 1. Let a be a right edge and assume it is regular. Writing a = g{c), let ip{N) = 
max{j : > c}. Then, almost surely, 

x^{N) —-> a, liminf (x,^(Ar)+i - a) > 0. (3.5) 

' ' N^oo N^oo ' ' 


Moreover, let {cn)n be the sequence associated with a as in Definition 3.1. Set 




\9'n{cn) 


1/3 


(3.6) 


aiv = 9n{(-n), cfn = 

so that aw —>■ a, cw —^ c, and ajsf —)• (2/y"(c))^/^ > 0 as N ^ oo. Then, for every s G M, 

Ji^P(^Ai^/^Crw(x<^(Ar) - ^n) ^ s) = lPAiry(ymax ^ s) • (3.7) 


Remark 3.8. Let us stress that the sequence (p{N) may be non-trivial even when considering 
the rightmost edge: As explained in Section 3.1, it is indeed possible that a certain amount 
of eigenvalues (possibly infinite) will converge outside of the limiting support. Thus, if we 
assume the rightmost edge is regular, then Theorem 1 states that there exists a maximal 
eigenvalue x,^(w) which actually converge to the rightmost edge and fluctuates according to 
the Tracy-Widom law. 

Let us comment on the history of this theorem. The Tracy-Widom fluctuations have been 
first obtained by Johansson [34] for the maximal eigenvalue when Hw is the identity. Baik, 
Ben Arous and Peche [7] then proved this still holds true when Sw is a finite rank perturbation 
of the identity, provided the perturbation is small enough. Assuming a condition which is 
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equivalent to the regularity condition (3.1) and that the maximal eigenvalue converges towards 
the rightmost edge, El Karoui [26] established the Tracy-Widom fluctuations for the maximal 
eigenvalue for general Sat’s assuming 7 ^ 1 , and Onatski [48] got rid of the last restriction. 
The statement on the existence of extremal eigenvalues converging to each regular right edge 
is [30, Theorem 2] and essentially relies on the exact separation results of Bai and Silverstein 
[5, 6 ]; the dehnition of the sequence ^p{N) indexing these extremal eigenvalues relies on these 
results. Finally, the Tracy-Widom fluctuations for the extremal eigenvalues associated to any 
right regular edge is [30, Theorem 3-(b)]. 

We now provide a similar statement for the left soft edges. 


Theorem 2. Let a he a left edge of the bulk. 7/ 7 > 1 and a is the leftmost edge of the bulk, 
set g:>{N) = n — N + 1. Otherwise, assume that a > 0 is regular, write a = g{c) and set 
ip{N) = min{Aj : < c}. Then, almost surely, 

x^{N) —-^ a, liminf (a - x^(at)_i) > 0. (3.8) 

^ ' N^OO N^OO ■r\ , , 

Moreover, let {cn)n be the sequence associated with a as in Definition 3.7. Set 


O-N = fl'iv(CAr), (Tat 



(3.9) 


so that Oat —>• a, cat —)• c, and (~ 2 / 5 "(c))^/^ > 0 as —>• 00 . Then, for every s G M, 


lim P 

N^OO 


(^N ^ ^ — IP^Airy (2/max ^ • 


(3.10) 


Prior to this result, which is a combination of Theorem 2 and Theorem 3-(a) from [30], 
the Tracy-Widom fluctuations for the smallest random eigenvalue when El at is the identity 
has been obtained by Borodin and Forrester [15]. 

Let us also mention that when v is the sum of two Dirac masses, a local uniform conver¬ 
gence to the Airy kernel (which is a weaker statement than the Tracy-Widom fluctuations) 
at every (regular) right and soft left edges follows from [41], see also [43, 44]. 

Finally, we state our last result, concerning the asymptotic independence of the Tracy- 
Widom fluctuations at a finite number of regular soft edges. For a more precise statement, 
we refer to [30, Theorem 4]. 

Theorem 3. Let be a finite collection of soft edges, and assume all these edges are 

regular. For each j G J, consider the rescaled eigenvalue N‘^^^crN,j{xip^(N) ~ ^N,j) associated 
with the soft edge a-j provided by (3.5)-(3.6) and (3.8)-(3.9). Then the family of random 
variables — aN,j)}j^j becomes asymptotically independent as N ^ 00 . 

The asymptotic independence has been previously established for the smallest and largest 
eigenvalues when El at is the identity by Basor, Chen and Zhang [9]. 

Remark 3.9. The results presented in this survey rely on the fact that the entries of Xat 
are complex Gaussian random variables, a key assumption in order to take advantage of the 
determinantal structure of the eigenvalues of the model under study. A recent work [37] 
by Knowles and Yin enables to transfer the results of Theorems 1, 2 and 3 to the case of 
complex, but not necessarily Gaussian, random variables. Indeed, by combining the local 
convergence to the limiting distribution established in [37] together with Theorems 1, 2 and 
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3, one obtains Tracy-Widom fluctuations and asymptotic independence in this more general 
setting, provided that the entries of matrix X^r fulfill some moment condition. Let us stress 
that the case of real Gaussian random variables (except the largest one covered in [40]), of 
important interest in statistical applications, remains open. 

We now turn to the hard edge and the Bessel point process. 


3.4 The Bessel point process at the hard edge 

The Bessel point process Tj^gsgi of parameter a G Z is the determinantal point process on 
M-i_ associated with the kernel 




^/y Ja{Vx)Ja{^) - Vx JaiVx)Jaiy/y) 


2 {x-y) 

where the Bessel function of the first kind Jq with parameter a is defined for x ^ 0 by 


(3.11) 




(-ir 


^n!r(n + a + l) V2/ 


n=0 


X\ 


2n 


(3.12) 


and satisfies the differential equation x^/"(x) + xf'{x) + (x^ — a^)/(x) = 0. 

The configurations (y^) generated by the Bessel point process a.s. have an infinite number 
of particles y* but have a smallest particle ymin- The law of ymin is characterized, for every 
s > 0, by 


^ [(ymin ^ s) — IP’BeLel((2^*) '®) — — det(I KgesseOi^fO.s)• (3.13) 


Bessel' 


When 0 = 0, this reduces to an exponential law of parameter 1, namely lPBgsggi(ymin ^ 
s) = e“®, as observed by Edelman [24]. In the general case, Tracy and Widom obtained the 
representation [52], 


(ymin ^ s) = exp ( / (logs - logx)g(x)^dx ) , 


p(“) 

^Bessel ^ 

where q is the solution of a differential equation which is reducible to a particular case of the 
Painleve V equation (involving a in its parameters) and boundary condition q(x) JaiVx) 
as X —>■ 0+. 

Recalling the Xj ’s are the eigenvalues of S^v and their distributional limit v has a compact 
support in (0,oo), we now provide our statement concerning the eigenvalues local behavior 
around the hard edge. 

Theorem 4. Assume that n = N + a with a G Z independent of N and set 

// . N 4^1 ^ 4 ,,,, , 8^1 


j=i J 


N^X) 

] = l 3 


Thus cjat —)• 4 J a ^z^(dA) > 0 and fjv —>■ 8 f X ^iy(dX) > 0 as N —>■ oo. 

Let Xmin be the smallest random eigenvalue o/Mat. Then, for every s > 0, we have 


Ji^P(^iV^CTArXmm ^ sj = PB^ggggi (ymin ^ s) ■ 


(3.14) 
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Furthermore, we have the expansion as N ^ oo, 



(3.15) 

When Sat is the identity, the convergence (3.14) has been established by Forrester [28]. 
As for the next order term (3.15), it has been obtained when Hat is the identity by Ferret 
and Schehr [49] and Bornemann [14], motivated by a question raised by Edelman, Guionnet 
and Peche in [25]. The statement (3.14) has been first obtained in [30], while the stronger 
statement (3.15) is [31, Theorem 6 ]. 


3.5 Application to condition numbers 


In this subsection, we study the fluctuations of the ratio 

3^ max 

Kn = - 

3Jmin 

of the largest to the smallest random eigenvalue of Mat. Notice that if re ^ A, then kn is 
the condition number of Mat while if re ^ A, then Kjsf is the condition number of Mat. The 
condition number is a central object of study in numerical linear algebra [47, 29]. Using our 
previous results, we can obtain an asymptotic description for kn- Let us emphasize that the 
leftmost edge a of the support of p is positive if and only if 7 / 1, see [31, Proposition 3]. 


Proposition 3.10. Assume 7 7 ^ 1. Denote by 0 the leftmost edge and by b the rightmost one. 
Assume that a, b are regular, Xmin —a and Xmax — >• b a.s. as A — )■ 00 (that a is regular and 
Xm\n —^ Cl a.s. is always true when 7 > Ij. Write a = g{c), b = g{d), consider the sequences 
(ctv) and (t)Af) associated with c and d respectively (see Definition 3.7) and set 


= gN{(-N) , 


Then, 



biv = 9N{bN) , 



kn 


a.s. 

-^ 

N—^00 


b 

a 


and 



V X bY 

-^ 7—I- 9 

Ai-s-oo oa cra^ 


T> 

where —> stands for the convergence in distribution and where X and Y are two independent 
random variables with the Tracy-Widom distribution. 


We now handle the case where 7 = 1 . 


Proposition 3.11. Assume n = N + a, where a € Z is independent of N, and moreover 
X ma x —)• b a.s. for some b > 0. Then, 


1 D 4b 

N^oo' ~x 


A ^re(dA) 


where P(A s) = Te^ssei(i/min ^ ■s) for every s > 0. 

Remark 3.12. Interestingly, in the square case where 7 = 1 , the fluctuations of the largest 
eigenvalue Xmax have no influence on the fluctuations oi kn as these are imposed by the 
limiting distribution of Xmin and the a.s. limit b of Xmax- 


Finally, we turn to the eigenvalues local behavior near a cusp point. 
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3.6 The Pearcey kernel at a cusp point 

Given any r G M, following [53] we introduce the Pearcey-like integral functions 

(^(x) = — / ^ r“ 

2^v^ /s 2z7r 

where the contour E consists in two rays going from ±e®’^/^oo to zero, and two rays going 
from zero to They satisfy the respective differential equations 

- T(p'{x) + x^p{x) = 0, (y) - 'Ttp'iy) - yV’(y) = 0. 

(t) 

The Pearcey point process Tpgai-cey determinantal point process associated with the 

Pearcey kernel 


K 


PeLey(®>y) = 


If”{x)'4){y) - Lp'{ x)'4j'{ y) + ip{x)Tp”{y) - T^p{x)Tp{y) 


x-y 


(3.16) 


This process has been first introduced by Brezin and Hikami [16, 17] when r = 0, and 
subsequent generalizations have been considered by Tracy and Widom [53]. 

The configurations (?/*) generated by the Pearcey point process are a.s. infinite and do 
not have a largest nor smallest particle. With this respect, the quantities of interest here are 
the gap probabilities of the Pearcey point process, defined for every 0 < s < t by 


ilrcey((yi) ^ [s,t] = 0 ) = det(/ - Kj,;^arcey)L 2 (.,*)• 


(3.17) 


Seen as a function of s, t and r, the log of the righthand side of (3.17) is know to satisfy a 
system of PDEs, see [53, 11, 1]. 

In [31, Theorem 5], we prove the following statement. 

Theorem 5. Let a = g{c) be a cusp point such that c G D, and assume it is regular. Let (cat) 
be the sequence associated with o as in Definition 3.7. Assume moreover the following decay 
assumption holds true: There exists k G M such that 


^ yiv(CAr) ^ 


N^oo 


(3.18) 


We set 


1/4 


<^N = y7v(Civ), 


<7AT = 


.aPicN), 


T = —K 




1/2 


(3.19) 


so that Oat —)• a and —)• {%/g'"> 0 as N ^ oo. Then, for every s > 0, we have 


^aN{xi - Oat)) n [-S, s] = 0 ) = ((yi) n [-s, s] = 0 ), 

where the Xi’s are the random eigenvalues of Mat . 


lim P 

N^OO 


(3.20) 


This result has been obtained by Mo when Sat has exactly two distinct eigenvalues [44] . 

As advocated in Section 3.2, the precise decay for ^ does influence the eigenval¬ 

ues local behavior near a cusp (see Proposition 3.6 and the discussion below). Our assump¬ 
tion (3.18) covers the general case where \^g'^{cN) —^ 0, and hence the limiting kernel is 
Kpgg^j,cey(a:, y) introduced by Brezin and Hikami, and the limiting regime where VNg'^{cN) 
has a limit as well. 
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Remark 3.13. (erosion of a valley) In the case where this limit k in (3.18) is positive, 
the deterministic equivalent measure /iv will not feature a cusp but rather a valley that will 
become deeper as —)■ oo, see the thin curve in Figure 5. The density of fiN will always be 
positive near the cusp and the condition 

5jv(C7v) ~ ^ 

should be thought of as a speed condition of the erosion of the valley. 

Remark 3.14. (moving cliffs) In the case where k < 0 in (3.18), 5^(cAr) is always negative 
for N large enough. In particular, there exists a small A'-neighborhood of Cjy whose image 
by is outside the support of /tat: There is a small hole in the support of but the 
two connected components move towards one another (moving cliffs), see the dotted curve in 
Figure 5. In this case, the condition 

ff7v(C7v) ~ ^ 

can also be interpreted as a speed condition at which the cliffs approach one another. 



Figure 5: Zoom of the density of /iv near the cusp point a. The thick curve is the density of 
^ in the framework of Figure 4. The thin curve (resp. the dotted curve) is the density of 
when VNg'i^{cN) > 0 (resp. VNg'^{cN) < 0). 

Remark 3.15. (slow decay) The slow decay setting where 

VNg'^{cN) -t Too 

is not covered by our results. In this case, we do not expect the Pearcey point process to arise 
anymore, and refer to Section 5 for further discussion. 
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4 Sketches of the proofs 

In this section, we provide an outline for the proofs of the results presented in Section 3. 


4.1 The random eigenvalues of Mat form a determinantal point process 

The key input on which all our proofs are based on, is that when the elements of X^r are 
complex Gaussian (Assumption 1), the configuration of the random eigenvalues Xj’s of Mat 
form a determinantal point process with an explicit kernel. More precisely Baik, Ben Arous 
and Peche provided in [7] a formula for that kernel, to which they give credit to Johansson. 
It is given by the following double complex integral 


^N{x,y) 


N 

(2i7r)2 



/ dy; 1 ^-Nx(z-q)+Ny{w-q) 

Je w- z 



(4.1) 


where the q G M is a free parameter (see [30, Remark 4.3]) and we recall the Aj’s are the eigen¬ 
values of Sat. The contours T and 0 are disjoint and closed, both oriented counterclockwise, 
such that r encloses all the A~^’s whereas 0 encloses the origin. 



Figure 6: The contours of integration 

Remark 4.1. The main ingredient to obtain this determinantal representation is the Harish- 
Chandra-Itzykson-Zuber integral formula, which allows to write a particular integral over the 
unitary group in terms of determinants, see [7, Section 2.1]. The analogue of this integral 
formula does not seem to exist for correlated Wishart matrices with real or quaternionic 
entries, and thus the determinantal structure seems only available in the complex setting. 

If we consider a random configuration of the form 

{N^(TN{Xi - Uat)) , 

then a change of variables yields it is a determinantal point process with kernel 

where Kat is as in (4.1). Hence, the study of the eigenvalues local behavior boils down to 
the asymptotic analysis as —)• oo of kernels of the form (4.2) with different choices for the 
scaling parameters /?, (Tat, uv- 
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4.2 Modes of convergence 

In order to prove the convergence (3.14) at the hard edge, it is enough to establish a local 
uniform convergence on M_|_ x M_|_ for the kernel (4.2) to the Bessel kernel y), after 

choosing appropriately the scaling parameters. Similarly, the local uniform convergence on 
MxM to the Pearcey kernel Kp’25,rcey(^’ n) yields the convergence (3.20) for the gap probabilities 
around a cusp. The convergences (3.7) and (3.10) to the Tracy-Widom law however require 
a stronger mode of convergence (such as the trace-class norm convergence, or the Hilbert- 
Schmidt norm plus trace convergence, for the associated operators acting on L^(s,oo), for 
every s G M; we refer to [30, Section 4.2] for further information). This essentially amounts 
to obtain a local uniform convergence on (s, -|-oo) x (s, +oo) plus tail estimates for K]y{x, y). 

From now, we shall disregard these convergence issues and provide heuristics on why the 
Airy kernel, the Bessel kernel and the Pearcey kernel should appear in different scaling limits. 

4.3 Towards the Airy kernel 

Here we provide an heuristic for the convergence to the Airy kernel. The gap to be filled in 
order to make this sketch of a proof mathematically rigorous can be found in [30] ; this heuristic 
may actually serve as a roadmap for the quite lengthy and technical proof we provided there. 

Since we are dealing with contours integrals, it is more convenient to use the following 
alternative representation for the Airy kernel (3.3), 

KAiry(x,y) = [ dz [ dn;—e—( 4 , 3 ) 

( 2 ^ 7 ^)^ w- z 

which is based on the contour integral formula for the Airy function (see e.g. the proof of 
[30, Lemma 4.15]). The contours H and E' are disjoint and unbounded contours, and E goes 
from to e~™/‘^oo whereas E' goes from to as shown on Figure 7. 



Figure 7: The paths of integration for the Airy kernel 


Consider the scaling parameters associated with a soft edge provided in 
thus focus on the right edge setting, but the situation for a left edge is similar. 


Theorem 1; we 
More precisely. 
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by using the formula (4.1) where we take q = cat, we investigate 


-Kjv( 


+ 


X 




AtVs 


dz 


I + 


1 


N'^/^aN 


(4.4) 


{2mYaN It 

where we introduced the map 


drc 
0 w- z 


^-N''^/3x(z-CM)/^N+N''-/^y{w-CN)/(^N ^NfN{z)-NfN{w) 


1 X ^ 

fN{z) = -aN{z - Cat) + log( 2 :) - log(l - Xjz). (4.5) 

J =1 

After performing the change of variables z cn + a^z/N^/^ and re i—)■ cat + the 

right-hand side of (4.4) becomes 

_ I _ / dz J) dte_ ^ _ ^-^^+y^+NfN{cN+<TN^:^)-NfN[tM+<rNj:^) gN 

( 2 i 7 r )2 Jvpfie) w- z 

where we set for convenience ‘PNiz) = N'^^^{z — cn)/(^n- 

Next, recalling gi^ was introduced in (3.2), the crucial observation that 

fN{z)=gN{z)-aN (4.7) 

allows to infer on the local behavior of /at around cat. More precisely, since by definition of 
the scaling parameters we have 

5Af(CAf) = <^N, 5jv(CAf) = 0, 9Ni'^N) -1 g"i'^) > 0, 

N^oo 

a Taylor expansion for /at around cat yields the approximation 

A^(/jv(civ+ (^jv^^) -/iv(civ)) ^ ^z^ , (4.8) 

where the constant 1/3 comes from the definition of ajsf, see (3.6). In conclusion, after plugging 
(4.8) into (4.6), we obtain the approximation 


Ar2/3 


(Tat 


Kat aN + 


X 


Ar2/3 


aN 


< flTV + 


iV2/3 


aN 


(2i7r) 


1 


dz f dm 
(/3jv(r) w — z 


—xz+yw+z^jS—w^fS 


(4.9) 


and we can almost read the Airy kernel (4.3), up to contour deformations. 

To frame the previous heuristic into a rigorous mathematical setting, a few technical points 
should be addressed, since of course the approximation (4.8) is only valid when \z\ is not too 
large and the contours appearing in the Airy kernel are unbounded. With this respect, the 
standard move is to split the contours T and 0 into different parts and then to deform each 
part in an appropriate way. 

In a neighborhood of cat, after simple transformations, one chooses T and 0 to match with 
the contours of the Airy kernel there, and then justify rigorously the approximation (4.9) after 
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restriction of z, w to that neighborhood. This can be done by quantifying the approximation 
(4.8) and then performing tedious but rather simple computations. 

Then, outside of this neighborhood, one proves that the remaining of the integrals don’t 
contribute in the large N limit. In the present setting of a general matrix XItv, this is the 
hard part of the proof. To do so, one establishes the existence of admissible deformations for 
the contours T and 0 so that they complete the Airy contours truncated on a neighborhood 
of Cat, and where the contribution coming from the term exp{Nfj^{z) — Nf]^{w)} brings 
an exponential decay on the remaining part. This can be done by looking for the so-called 
steepest descent/ascent contours (i.e. contours on which Re /tv is decreasing/increasing), and 
this was the strategy used by Baik, Ben Arous, Peche [7] and El Karoui [26] when dealing 
with the rightmost edge. When considering any right or left soft edge, following this strategy 
requires to consider many sub-cases and to perform again most of the computations in several 
case. In [30], we instead developed a unified (abstract) method to provide the existence of 
appropriate contours by mean of the maximum principle for subharmonic functions. 

For the reader interested in having a look at the proofs of [30] , let us mention it turns out 
it is more convenient to work at a scale where the contours T, 0 live in a bounded domain, and 
this is the reason why we did not performed there the changes of variables z i—)• CN + a^z/N^^^ 
and tc I—)■ Ctv + ai\fw/N^/^ as we did in the present heuristic. 


4.4 Towards the Pearcey kernel 

Now we turn to the heuristics for Pearcey kernel and refer to [31] for a rigorous proof. The 
setting is essentially the same as in the Airy case, except that now ctv is a simple zero for 
instead of g'j^, and g'j^icjy) —)• 0, which entails a different behavior for the map /tv near ctv- 
We start with the alternative representation for the Pearcey kernel (3.16), 


Pearcey 


{x,y) 


1 

(2i7r)2 


dz 


dw 


2 4 2 4 


w — z 


(4.10) 


where the contour H is disjoint to the imaginary axis, and has two components. The first 

part goes from e^'^/^oo to whereas the other part goes from to 

See [53] for a proof (and also [16] when r = 0). Notice also the symmetry ^Peaiceyi^^y) ~ 

Kj/gg^j,cey(— X, —y) which follows from the change of variables z,w i— —z,—w. 

Consider the scaling parameters associated with a regular cusp point provided in Theo¬ 
rem 5. By using the formula (4.1) where we choose q = ctv, we now consider 


iVl/4 r ^ r 1 _Ari/42,P^Ov)+TVi/4^hf!^M) 

= 7 -TTT- (b dz (t dw -e ‘’N <^N 

{2nTyaN Jr Je w- z 




(4.11) 


where the map /tv is the same as in (4.5). After the change of variables z i—>• ctv + ajsizlN^^^ 
and tc I—)■ Ctv + cr^w, the right-hand side of (4.11) reads 


1 

(2i7r)2 


dz ® dtc 
‘Piv(r) 4(pjv(©) 


1 ^-xz+yw+Nfpf (civ+o-jv )-^fN ) 

W — Z 


(4.12) 
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where we introduced for convenience ^n{z) = — C]sf)/o'N- 

In this setting, the definition of the scaling parameters yields 

ffiv(CAr) = aw, ViV5^(cAr) —- > K, ^^(cAr) = 0, 5w(CAf) T;-^ /'(c) > 0. 

N—>-oo N—^oo 

Recalling the identity (4.7) and the definition (3.19) of r, a Taylor expansion around cw then 
yields the approximation 

A^(/7v(cw + CTw^^) -/7v(cw)) ^ + 

where the constant 1/4 comes from the definition of ajy, see (3.19). Thus, by plugging (4.13) 
into (4.12), we obtain the approximation 


Ar3/4 


(Tat 


Kat 1 aw + 


X 


Ar3/4 


fJw 


I aw + 




(2z7r)^ 


1 


dz f dw 
v^ivT) JfNie) ^ ^ 


2 4 2 4 

-xz-^ + ^+yw+^ ^ 


(4.14) 


and we can almost see the Pearcey kernel (4.10), up to contour deformations. 

As to make this approximation rigorous, the method is the same as for the Airy kernel. Let 
us mention that the abstract argument we mention previously for the existence of appropriate 
contour deformations also applies in this setting. 


4.5 Towards the Bessel kernel 


Finally, we provide heuristics for the appearance of the Bessel kernel and refer to [31] for a 
rigorous proof. The main input here is, according to Section 2, the critical point c associated 
with the hard edge is now located at infinity. 

The first step is to write the Bessel kernel (3.11) as the double contour integral. 


K 


(q) 

Bessel 


{x,y) 


1 

(2i7r)^ 


f I ^ f z\°‘ 

x) J\z\ = r Z J\u,\=R w z-w\w) 


2 

4 


y. 

w 


(4.15) 


where 0 < r < i?, and which is provided in [30, Lemma 6.2]. The contours of integration 
are circles oriented counterclockwise. Let us stress this formula is only available when a G Z, 
since otherwise the term (z/w)°‘ in the integrand would not make sense on the whole of the 
integration contours. 

Setting fJw as in Theorem 4 and using the formula (4.1) where we choose q = 0, we now 
consider 


A'^dw ( A^dw ’ A^dw) 


1 

(2f7r)2Adw 




N zx 

e 



1=1 


w - A-^ 


(4.16) 


Having in mind the critical point is located at infinity, we perform the change of variables 
2:1—)• Nctn/z and w 1—)■ Na^/w, so that the right-hand side of (4.16) reads 


1 

(2f7r)2 


/ 

4 ‘/’ivfr) 


2 


/ 


dtc 

w 


1 


z — w 



1-NGn(z)+NGm{w) 

1 


(4.17) 
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where we introduced the maps 


and (Pn{z) = Na^/z. We emphasize that, during the previous step, we used that n = N + a 
and witnessed a cancellation leading to the term {z/w)°‘, which does not depend on N. 

Now, a Taylor expansion of Gw around zero yields the approximation 

iv(Gw(^)-Gw(0)) c. . (4.18) 

Thus, by plugging (4.18) into (4.17), we obtain the approximation 
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Kw 


( ^ y \ 

\N^aN' N‘^(Tn) 

1 f dz f dw 1 


( 2 i 7 r )2 z J^^( 0 ) w z-w 



y. 

w 


T, (4.19) 


and we can almost see the Bessel kernel (4.15), up to contour deformations and the pref¬ 
actor (y/x)“/^. Finally, in order to deal with that prefactor, one considers the operator E 
of multiplication by x“/2 acting on L2(0 ,s), and then use that replacing the Bessel kernel 
K^gsei(®)y) by the kernel of the operator which is (4.15) without the prefactor 

(y/x)"/2^ leaves the Fredholm determinant det(/ — Ki^sgei)i^(o,s) invariant. 

To make this heuristic rigorous, the method is far less demanding than in the setting 
of the Airy or the Pearcy kernel. Indeed, in the present setting one can legally deform the 
contours F and 0 in such a way that (^w(r) and (/?Ar(0) match the integration contours for 
the Bessel kernel (4.15). After making this deformation, a simple Taylor expansion of the 
map Gw around zero will be enough to establish the convergence towards the Bessel kernel 
and therefore to obtain (3.14). 

By pushing the Taylor expansion (4.18) one step further, one can also obtain the more 
accurate estimate (3.15), provided with an identity involving the resolvent of the Bessel kernel 
established by Tracy and Widom [52]. We refer the reader to [31] for further information. 


5 Open questions 

The results presented here naturally entail a number of open questions that we list below. 

1. At the edge of the definition domain and exotic local behaviors. The results 
on the eigenvalues local behavior presented in this survey only concern edges or cusp 
points a which read a = ^(c) with c € D. If we focus on the rightmost edge for the 
sake of simplicity, then Proposition 2.1 states this edge may actually belong to g{dD) 
(notice that this cannot happen if the limiting spectral distribution u of Sw is a finite 
combination of Dirac measures). In this case, the square root behavior of the density 
around this edge is not guaranteed anymore, and the laws describing the fluctuations of 
the eigenvalues near such an edge seem completely unknown and a priori different from 
Tracy-Widom distribution. We believe the fluctuations will actually depend on z^, and 
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hence lie outside of the random matrix universality class. Quite interestingly, the same 
phenomenon arise in the study of the additive deformation of a GUE random matrix 
[20] and random Gelfand-Tsetlin patterns [22, 23]. 

2. Alternative regime at a cusp point I. In the context of Theorem 5, our speed 
assumption (3.18) does not cover the following case 

^9n{^n) —-^ +00 . 

N^oo 

This condition corresponds to the situation where the density of the deterministic equiv¬ 
alent /Utv is positive in a neighborhood of cv- It essentially states that the bulk of /iv 
will degenerate into a cusp around g{c) quite slowly and we do not expect to witness 
Pearcey-like fluctuations around S'Ar(cAf) anymore. We believe instead that the sine 
kernel will arise at the scale \/N/g'j^{cN), which strictly lies in between and 

3. Alternative regime at a cusp point II. Another case that is not covered by our 
assumption (3.18) is when 

^ 5jv(cv) —— -)> —oo . 

N^oo 

In this case, cv lies outside the support oi fiN and g'j^ has two distinct real zeroes 
near cv, say cv,! and CAr,2- Hence, for N sufficiently large, gN{(-N,i) and 5Ar(cv,2) both 
correspond to edges of the support of /xjv which both converge towards the cusp point 
g{c). The previous condition entails that such a convergence will happen at a quite 
slow rate. In this case, we do not expect to observe the Pearcey kernel around ^^(cv) 
either, because of the absence of particles, but we believe a local analysis around the 
edge gN{(-N,i) or 5Ar(cv,2) niay uncover the Airy kernel at an intermediary scale. 

4. Study of the fluctuations at the hard edge in more general cases. The hard 
edge fluctuations were described here when n = N + a with a fixed a G Z , but the hard 
edge is always present as soon as re/A —)• 1. Thus it would be of interest to describe the 
hard edge fluctuations in more general situations, for example when a = a (re) —>• -|-cx) 
so that re/A —)• 1. In this case one would expect Tracy-Widom fluctuations near the 
leftmost edge of fijy, the latter being positive but converging to zero as A —)■ oo. 

5. Non-Gaussian entries. 

All the fluctuations results presented here rely on the fact that the entries of matrix 
Xtv are complex Gaussian. It is however of major interest, for applications and for 
the general theory as well, to study the universality of such results for non-Gaussian 
complex random variables. As explained in Remark 3.9, the Tracy-Widom fluctuations 
for the extremal eigenvalues associated with any regular soft edges are now established 
in this general setting (under some moment conditions for the entries), by combining 
theorems 1 and 2 with Knowles and Yin’s recent preprint [37] (see also [8]). However, 
natural related questions remain open: Would it be possible to describe the fluctuations 
at 

(a) the hard edge, for general complex entries? 

(b) a (regular) cusp point, for general complex entries? 
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Another universality class of interest is the case where the matrix Xjv has real entries. 
In this case, the techniques based on the determinantal structure of the eigenvalues 
are no longer available. Lee and Schnelli [40] recently succeeded to establish GOE 
Tracy-Widom fluctuations of the largest eigenvalue when the entries are real Gaussian 
or simply real (with subexponential decay), under the assumption that the covariance 
matrix Xtv is diagonal. The techniques developed by Knowles and Yin [37] enable to 
relax the diagonal assumption for the covariance matrix SItv. A number of questions 
remain open: Would it be possible to describe the fluctuations at 

(c) any (regular) soft edge, when the entries are real Gaussian? 

(d) the hard edge, when the entries are real (Gaussian or not)? 

(e) a (regular) cusp point, when the entries are real (Gaussian or not)? 
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